O(N) Iterative and O(NlogN) Fast Direct Volume

Linear complexity iterative and log-linear complexity direct solvers are developed for the volume integral equation (VIE) based general large-scale electrodynamic analysis. The dense VIE system matrix is first represented by a new clusterbased multilevel low-rank representation. In this representation, all the admissible blocks associated with a single cluster are grouped together and represented by a single low-rank block, whose rank is minimized based on prescribed accuracy. From such an initial representation, an efficient algorithm is developed to generate a minimal-rank H-matrix representation. This representation facilitates faster computation, and ensures the same minimal rank’s growth rate with electrical size as evaluated from singular value decomposition. Taking into account the rank’s growth with electrical size, we develop linear-complexity H-matrix-based storage and matrix-vector multiplication, and thereby an O(N) iterative VIE solver regardless of electrical size. Moreover, we develop an O(NlogN) matrix inversion, and hence a fast O(NlogN) direct VIE solver for large-scale electrodynamic analysis. Both theoretical analysis and numerical simulations of large-scale 1-, 2and 3-D structures on a singlecore CPU, resulting in millions of unknowns, have demonstrated the low complexity and superior performance of the proposed VIE electrodynamic solvers.